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The basic idea of a microscopic understanding of Thermodynamics is to derive its 
main features from a microscopic probability distribution. In such a vein, we investigate 
the thermal statistics of quasi-probabilities’s semi-classical analogs in phase space for 
the important case of quadratic Hamiltonians, focusing attention in the three more 
important instances, i.e., those of Wigner, P-, and Husimi distributions. Introduction of 
an effective temperature permits one to obtain a unified thermodynamic description that 
encompasses and unifies the three different quasi-probability distributions. This unified 
description turns out to be classical. 

Keywords: Semiclassical physics; P—function; Information quantifiers; thermal proper¬ 
ties 


1. Introduction 

A quasi-probability (QP) distribution is a mathematical construction resembling 
a probability distribution. It does not necessarily fulfill, though, some of the Kol¬ 
mogorov axioms for probabilities QPs display general features of ordinary proba¬ 
bilities and yield expectation values with respect to the weights of the distribution. 
However, they disobey the third probability postulate 12, in the sense that regions 
integrated under them do not represent probabilities of mutually exclusive states. 
Some quasi-probability distributions exhibit zones of negative probability density. 
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QPs often arise in trying to study quantum mechanics a phase space representation. 
This is frequently done in quantum optics, time-frequency analysis, etc. 

The dynamics of a quantum system is determined by a master equation. We 
speak of an equation of motion for the density operator (p), defined with respect to 
a complete orthonormal basis. One can show that the density can always be written 
in a diagonal manner, with re spec t to an overcomplete basis If this basis is that 
of coherent states |a) one has^^Ml 

^ = y (1) 

where d^a/7r = dxdp/2Trh, with x and p variables of the phase space. The sys¬ 
tem evolves as prescribed by the evolution of the quasi-probability distribution 
function. Coherent states, eigenstates of the annihilation operator d, serve as the 
overcomplete basis in such a build-up ^2121 

There exists a family of different representations, each connected to a different 
ordering of the creation and destruction operators a and d^. Historically, the first of 
these is the Wigner quasi-probability distribution IF® related to symmetric opera¬ 
tor ordering. In quantum optics the particle number operator is naturally expressed 
in normal order and, in the pertinent scenario, the associated representation of the 
phase space distribution is the Glauber-Sudarshan P one In addition to IF and 
P, one may find many other quasi-probability distributions emerging in alternative 
representations of the pha se space distribution®. A quite popular representation is 
the Husimi Q one I ^ I ^ I S ISI^ used when operators are cast in anti-normal order. 

In this paper we wish to apply semiclassical information theory tools associated 
to these P, Q, and IF representations (for quadratic Hamiltonians) in order to 
deseribe the concomitant thermal semiclassical features. We specialize things to the 
three /—functions associated to a Harmonic Oscillator (HO) of angular frequency 
w. In such a scenario the three functions -that we name for sake of convenience fp, 
fg, and fw- are simple Gaussians and the treatment becomes entirely analytical, 
a very convenient feature. The HO is a really important system that yields insights 
usually having a wide impact. Thus, the HO constitutes much more than a mere 
simple example. Nowadays, it is of particular i nterest for the dynamics of bosonic 
or fermionic atoms contained in magnetic traps [IHIlIIISl ^g for any system 

that exhibits an equidistant level spacing in the vicinity of the ground state, like 
nuclei or Luttinger liquids. 

In this paper we are interested in thermal states regarded as semi-classical dis¬ 
tributions in phase space -analogs of the quantum quasi-probabilistic distributions. 
We will try to ascertain what physical Thermodynamics’ f eature s are they able to 
describe at the semi-classical level. These distributions are^^^HH 




^huj 


, 7P = e' 


I (P — function), 


( 2 ) 
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Iq = IQ 7 Q = 1 - (Q - function), (3) 

fw = Iw = 2 tanh(/l/kt;/2) (VF — function), (4) 

with /3 = l/ksT, ks the Boltzmann constant, and T the temperature. As stated 
above, these distributions will be used in the next section as semiclassical statisti¬ 
cal weight functions. Since ours is NOT a quantum approach, the ordering of the 
HO-creation and destruction operators a and plays no role whatsoever in our 
manipulations below. 

This paper is organized as follows: section [prefers to information quantifiers, in 
a phase space representation, for Gaussian distributions. In Section [3] we focus at¬ 
tention upon thermodynamic relations and we express them in terms of an effective 
temperature. Finally, some conclusions are drawn in Section |4l 


2. Semi-classical information quantifiers 

• The first step in our development is to calculate entropic quantifiers for 
these Gaussian distributions. In order to simplify the notation we will con¬ 
sider a general normalized gaussian distribution in phase space 

/(a) = 7e-^l“l', (5) 

whose normalized variance is I /7 and 7 taking values 7 p, 7 q and "fw given 
by Eqs. (ED, ED, and ED, respectively. 

The logarithmic Boltzmann’s information measure for the the probability 
distribution ED is 

S = -J ^ f{a) \nf (a) = l-\n'y, ( 6 ) 

so that it acquires the particular values 


5'p = I - In - I) , 

(7) 

5q = I - In (I - e-d^) , 

( 8 ) 

Sw = I — ln(2tanh(/3fia;/2)), 

(9) 


for, respectively, the distributions fp, /g, and fw- 
• Next, we focus attention on the information quantifier known as Fisher’s 
information measure. We specialize it for families of shift-inva riant d istri- 
butions, that do not change shape under translations. One has^^^ES! 


1 = 


dxf{x) 


/ d\nf{x)\ 

V ) 


I 


I 


2 


( 10 ) 
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and, in phase space, it adopts the appearance 


[m 




/ d\nf{a) \ 
\ d\a\ ) 


( 11 ) 


such that considering f{a) given by Eq. ® we get / = 7 , whose specific 
values are yp, yg, yw for the three functions fp, /g, and fw, respectively. 

• The statistical complexity is a functional C[P] that can primarily be viewed 
as a quantity that characterizes the probability distribution P. It quant ifies 
not only randomness but also the presence of correlational structures 
The opposite extremes of perfect order and maximal randomness possess 
no structure to speak of. In between these two special instances, a wide 
range of possible degrees of physical structure exist, degrees that should 
be reflected in the features of the underlying probability distribution. The 
statistical complexity C, according to Lopez-Ruiz, Mancini, and Calvet^^, 
is a suitable product of two quantifiers, such that C becomes minimal at the 
extreme situations of perfect order or total randomness. Instead of using 
the prescription of^^, but without violating its spirit, we will take one 
of these two quantifiers to be Fisher’s measure and the other an entr qpic 
form, since it is well known that the two behave in opposite manner ( 1 ^. 
Thus: 


C = 5/ = y(I-lny), (12) 

that vanishes for perfect order or total randomness. For each particular 
case, we explicitly have 


Cp = - 1) [1 - In - l)] , (13) 

Cg = (I - [I - In (I - e-^^)] , (14) 

Cw = (2 tanh(,5?ku/2)) [I — In (2 tanh(,5?ku/2))], (15) 

for, respectively, the distributions fp, /g, and fw- The maximum of the 
statistical complexity occurs when y = I and, the associated temperature 
values are 

' — I = I T = huj/ks In 2 for the /p—function, 

< I — = I => T = 0 for the /g—function, 

, 2 tanh(/3/ia;/2) = I => T = /ka/2/cp arctan(I/2) for the /w—function. 
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3. Thermodynamic relations 

We start this section considering the semi-classical Hamiltonian of the harmonic 
oscillator that reads 

= huj\a\'^ = fuv , (17) 

where x and p are phase space variables and cr^ = h/2mui and = hmui/2^^. Let 
us further define the semiclassical expectation value of the function A{x,p) as 

{A)f=J^f{a)A{x,p), (18) 

indicating that /(o) is the statistical weight function. Thus, the mean energy of 
the hamiltonian 77(x,p) is written in the fashion 

U* = hw [ = (19) 

7 TT 7 

where 7 takes the respective values 7p, 7 q, and explained in Introduction. 
Additionally, the thermodynamic entropy S is 

S =fcB(l-ln7), (20) 


where we have added the Boltzmann constant fcp. The mean energy can be viewed 
as a function of the thermodynamic entropy S , in the following fashion. Combining 
dm with the thermodynamic entropy (1201) we get the associated, fundamental 
equation U* = U*{S ) 


U*{s') = fiwe^'TB-fi 

(21) 

and 


7 = ei-s'AB. 

(22) 

Thus, the differential of U* becomes 



(23) 


where we have considered the volume V to be constant. Thus, after effecting the 
pertinent replacements we find 


dU* = f^dS', (24) 

kbj 

which suggests introducing an effective temperature Te//- Using Teff we obtain a 
unified picture that encompasses the three distributions fp, fg, and fw, in a single 
thermodynamic description. We have 


dU* \ _ 
dS' ) y fcsy ’ 


I 


I 


Teff = 


(25) 
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such that 

dU*=TeffdS'. (26) 

Note that in the three instances, Fe// = cxd for T = oo. However, if T = 0, 
Teff =0 only in the /p-case. It equals 1/2 in the Wigner instance and equals 1 in 
the Husimi case, as depicted in the accompanying figure. 



Fig. 1. Effective temperature Teff versus temperature T in 


From dm) and (j26|l we can rewrite the mean energy in terms of effective tem¬ 
perature. 


U*=kBT,ff, 


(27) 


that corresponds to the classical mean energy of a harmonic oscillator of temper¬ 
ature Fe//, with /cpFe///2 contributions for each of the two pertinent degrees of 
freedom. Similarly, the thermodynamic entropy is recast as 


^ = 1 + In 
Kb 


fcg Teff 

Hlo 


and the Helmholtz free energy is given by 

A* = U*-TeffS'=kBT,ff In 
The effective specific heat is defined as 

Cf = 


huj 


ks Teff 


dU* 


dT, 


eff/v 


(28) 


(29) 


(30) 


that using (EZD becomes 
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which is precisely the specific heat for the classical harmonic oscillator which is 
independent of the temperature. This becomes the Dulong and Petit’s rule at the 
classical limit. In view of ((26ll and ((2^ the analog partition function Z* is given by 


z* = 

(32) 

7 

and, according to Eqs. (jH), (IT^. and (1^ we find 


s' =\nZ*+P*U*, 

(33) 

with 


/3* — ^ 

(34) 

^ ksTeff hhj 


Thus, one reobtains all the thermal results pertaining to a classical HO at the tempera¬ 
ture Teff. Note that the whole thermal description becomes now of a classical character. 
All the quantum effects pertaining to the probability distributions are contained in the 
relationship (24) between Teff and T. 


It is interesting to look at the statistical complexity C in order to see at what ef¬ 
fective temperature the possible correlational structures carried by our probabilistic 
distributions are stronger. Expressed in terms of Te//, C becomes 


C' = IS = 


huj 

Teff 


1 -b In 


^ ksTeff '^ 


(35) 


Keeping in mind Te/z’s definition, it is easy to see that the maximum for the 
complexity C is attained when Teff = huj/kB- This implies, according to Eq. (1251) 
that the maximum of the Fisher measure es Imax = 1. At the complexity-peak, 
thermodynamic quantities take the values 


^max 

(36) 

^max ? 

(37) 

^max 1 

(38) 

a remarkable simplicity! 


4. Conclusions 


We have investigated here the thermal statistics of quasi-probabilities-analogs /(a) 
in phase space for the important case of quadratic Hamiltonians, focusing atten- 


tion on the three more important instances, i.e., those of Wigner, P-, and Husimi 
distributions. 
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• Introduction of an effective temperature permits one to obtain a unified 
thermodynamic description that encompasses the three different quasi¬ 
probability distributions. This unified description turns out to be classical. 

• The above entails that all possible “quasi-quantum” effects have to be 

contained in the relationship between T^ff and T. Note, for instance, that 
the minimal energy is not zero {one of these effects) but, in the /w—case, 
Tjjy = haj/2kB, implying a minimum energy = fujj/2. Additionally, 

the Husimi-Tg'j^ = huj/kB reflects the well known fact that the Husimi 
distribution “smoothes” the Wigner one over a phase-space area = h. 

• The basic idea of a microscopic understanding of thermodynamics is to 
derive its main features from a microscopic probability distribution. We 
have done just this using as distributions quasi-probability ones. 
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